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We study a binary Lennard-Jones system below the glass transition with molecular dynamics 
simulations. To investigate the dynamics we focus on events ("jumps") where a particle escapes 
the cage formed by its neighbors. Using single particle trajectories we define a jump by comparing 
for each particle its fluctuations with its changes in average position. We find two kinds of jumps: 
"reversible jumps," where a particle jumps back and forth between two or more average positions, 
and "irreversible jumps," where a particle does not return to any of its former average positions, i.e. 
successfully escapes its cage. For all investigated temperatures both kinds of particles jump and both 
irreversible and reversible jumps occur. With increasing temperature relaxation is enhanced by an 
increasing number of jumps, and growing jump lengths in position and potential energy. However, 
the waiting time between two successive jumps is independent of temperature. This temperature 
independence might be due to aging, which is present in our system. The times of jump duration 
are at all temperatures significantly shorter than the waiting times. The ratio of irreversible to 
reversible jumps is also increasing with increasing temperature, which we interpret as a consequence 
of the increased likelihood of changes in the cages, i.e. a blocking of the "entrance" back into the 
previous cage. In accordance with this interpretation, the fluctuations both in position and energy 
are increasing with increasing temperature. A comparison of the fluctuations of jumping particles 
and non-jumping particles indicates that jumping particles are more mobile even when not jumping. 
The jumps in energy normalized by their fluctuations are decreasing with increasing temperature, 
which is consistent with relaxation being increasingly driven by thermal fluctuations. In accordance 
with subdiffusive behavior are the distributions of waiting times and jump lengths in position. 

PACS numbers: 02.70.Ns, 05.20.-y, 61.20.Lc, 61.43.Fs, 64.70-Pf 
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I. INTRODUCTION 

If a liquid is cooled and crystallization is avoided, one 
obtains a Supercooled liquid. Upon further cooling the 
system falls out of equilibrium and results in a glass. 
During the transition from liquid to supercooled liquid 
to glass the thermodynamic properties change and even 
more drastic changes occur in the dynamics p|. The 
viscosity increases by many orders of magnitude upon 
cooling and the mean square displacement (MSD) as a 
function of time develops a plateau at intermediate times. 
The lower the temperature the longer the waiting time 
within the plateau until a second rise in the mean square 
displacement occurs 2]. One common explanation for 
the plateau is that while at high temperatures one has 
(at late enough times) normal diffusion, at lower temper- 
atures particles are caged in, i.e. trapped by their neigh- 
bors, and spend longer time within this cage with de- 
creasing temperatures. The second increase in the MSD 
indicates that after long enough waiting time the parti- 
cles manage to escape their cage. This escape out of the 
cage ("jump") is the focus of this paper. 

To set the work of this paper in context, we review 
briefly previous studies on the dynamics of supercooled 
liquids and glasses. We restrict ourselves mostly to sim- 
ple glass formers, excluding major work about strong 
glass formers, ionic systems, polymers and crystals. 



Central quantities in both experiments and com- 
puter simulations arc the viscosity, diffusion constant 
and MSD. The MSD, which is an average over single par- 
ticle i — 1, . . . ,N displacements. 
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where r^(t) is the position after steepest decent 

show jumps when viewed with fine enough time resolu- 
tion, and when studied at low enoiigh temperatures and 
for single configurations |llllEISIilElSll3- De- 
tailed studies of these jumps indicate collective motion 
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One learns about the interplay between regular diffu- 
sion and hopping motion via four-point correlation func- 
tions [Til Il2l ITsL 01 and the van Hove correlation 
function G{rA) = Gs{r,t) + Gd{r,t) [HIIllj?, 18, 1^ 
113, El m IM il where Gs{r,t) is the self part and 
Gdir, t) is the distinct part. The hopping shows up in an 
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additional peak in Gg (r, t) and an increase in the ampli- 
tude at r ^ of G^{r, t). 

Recently attention has been drawn to the non- 
Gaussian tail of G{r,t) [2 1 [ M | 27[ [2! and the non- 
Gaussian parameter a2 jlllllMlallsS IM El 1^ 
Isl Is^l which is the second coefficient of the cumulant 
series of Fs {q, t) , the Fourier transform of Gs {r, t) . They 
are signatures of non-exponential behavior, which might 
be either due to homogeneous complex dynamics and/or 
due to spatially heterogeneous dynamics 
been studied both with experiments 
in simulations of two-dimensional 
three-dimensional systems 0, |^ 

E El El El H El Studies of the most mobile 
particles reveal that they move increasingly more collec- 
tively with decreasing temperature. One distinct feature 
of this collective dynamics is that some particles are mov- 
ing along a string H i, H H El El El ■ 

Similar motion has been found via the study of normal 
modes which are commonly used to study the dynamics 
at low temperatures 0, El El El El El i^i] . 

The dynamics of glasses out of equilibrium, i.e. sys- 
tems which have been quenched from high to low tem- 
perature, displays additional complexity, since the sys- 
tem might "age" ^J. This means that the dynamics 
depends on the waitin g ti me between the quench and the 
measurement |5ll l52l l53j| which results in the violation 
of the fluctuation dissipation theorem [sl 0| . 

Another fruitful approach to gain insight into the dy- 
namics of supercooled liquids and glasses has been to in- 
vestigate the energy landscape of the inherent structure, 
that means of the instantaneous configurations which 
have been quenched to their local potential energy min- 
imum [11 iHl m m m m At the same temper- 
ature when the system starts slowing down drastically, 
the potential energy of the inherent structure undergoes 
a qualitative change [s^ . The long time behavior of this 
potential energy shows jumps between metabasins, where 
the latter are you ps of strongly correlated local min- 
ima [5l[5l[58', 59, 60]. The mean average waiting time 
(Atwait) between these jumps turns out to be dominated 
by the long times. (At^ait) together with the diffusion 
constant allow an estimate of the cage size [sl] . 

Direct studies of the cage have been done via three- 
time correlations 0|, velocity- velocity correlations [sl 
in [^ I via the cage correlation function "6?, , and 
a passage time before a particle escapes ,65j . The results 
of these studies tell us about cage properties such as the 
cage size and waiting time within a cage (or jump rate). 

Hand in hand with the experiments and simulations 
are the theoretical models for supercooled liquids and 
glasses. One of the very successful theories is the mode 
coupling theory of the glass transition (MCT) 66] , which 
describes the dynamics of supercooled liquids via non- 
linear dynamics of coupled density modes and makes pre- 
dictions for quantities such as D, F(q, t), and the suscep- 
tibility The extended version of MCT includes 
hopping processes via the coupling to current fluctuations 



[67|. The comparison of experiment and this theory for 
x"(w) [il] and a2 Q9\ shows very good agreement. 

The MGT is a theory for the glass transition for tem- 
peratures above the glass transition. Below the glass 
transition there exists no equivalent of a microscopic the- 
ory such as MGT. Some phenomenological theories which 
are more general models for solids build in hopping either 
indirectly as in the free volume model 70] or directly in 
hopping models [^ [71 [t^ • Direct hopping models are 
based on the assumption that hopping is the main process 
which explains the dynamics. The starting point in these 
models is usually a random walk and the anomalous dif- 
fusion is incorporated via a distribution of waiting times 
(or jump rates). The jumping entities are either defects 
fflf , uncharged carriers ^ or (in the case of predictions 
for ac and dc conductivity) charged carriers [73. 

The goal of this paper is to obtain via single particle 
trajectories direct information about jump processes such 
as jump sizes and waiting times between jumps. Specific 
examples of the dynamics of single particle jumps , in the 
form of a plot of one component of (t) [l^ [13, [H UM 
or in the form of a two- or three-dimensiona,l picture of 
the particle trajectory [1 12I I29I Isol Isil liollSTlleH . are 
very helpful to get an idea of some qualitative features of 
jumps, such as the detailed geometry of jump processes 
[H[H[H[ZJ|' In this paper, however, we go beyond sin- 
gle examples by defining a systematic search algorithm 
for jump processes. For the case of single particle tra- 
jectories a similar approach has been used in the work 
\VA [19I [SJ, [43, where the jump criterion is a minimum 
hopping distance, which is the same for all particles. In 
contrast, we use a relative criterion, where for each par- 
ticle its size of fluctuations is compared to its jump size. 
We choose this relative criterion to be able to identify 
jumps of particles of different sizes and neighborhood. 
That means the criterion is adjusted to the cage size of 
each individual particle. As criterion for the occurrence 
of a jump we use the positions of particles instead of 
their energy. While both approaches are fruitful, we be- 
lieve that many theoretical models and interpretations 
of simulations and experiments are based on an intuitive 
picture of the particle motion in real space. Similarly 
our goal is to mimic a careful observer of each particles 
motion in our system. We therefore use single particle 
trajectories, instead of a quantity which is an av erag e 
over all particles (as in the work [1 H H IllTlHlgllTof) . 
For the case of single particle trajectories the distinction 
between vibrational and hopping motion turns out to be 
clearer by taking time averages than by using trajectories 
of the inherent structure. 

In the following we define our model, and give details 
about the simulation (Sec. Our precise definition of 
a jump is given in Sec. IIIII We find two types of jumps 
(irreversible and reversible) . The latter are distinguished 
for the rest of the paper (in distinction from flTl Il9l [SJ, 
El .) In Sec. II VI we count as a function of temperature the 
number of jumping particles, and in Sec. the number 
of visited different cages. In Sec. IVII we investigate the 
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times during a jump and between successive jumps. The 
jump size both in position and in potential energy are 
presented in Sees. IVIll and Ivml In Sec.Hxl we conclude 
with a summary of our results, a comparison with the 
results of previous work and with our resulting picture of 
jmnp processes. We finish with open questions suggesting 
future work. 



II. MODEL AND SIMULATION 

We use a binary Lennard-Jones (LJ) mixture of 800 A 
and 200 B particles with the same mass. The interaction 
potential for particles i and j at positions and Vj and 
of type a,/3 e {A,B} is 

- 4.,., ((f^)" (5) 

where r = [r^ - r^l and eaa = 1-0, eab = 1-5, ebb = 
0.5,crAA = 1-0, (Tab — 0-8 and (Tbb = 0.88. We trun- 
cate and shift the pot ential at r = 2.5 ctq./?. From previ- 
ous investigations |75l it is known that this system 
is not prone to crystallization and dcmixing. In the fol- 
lowing we will use reduced units where the unit of length 
is (TAA, the unit of energy is eAA and the unit of time is 
V^mcr^A /(48eAA)- 

We carry out molecular dynamics (MD) simulations 
using the velocity Verlet algorithm with a time step of 
0.02. The volume is kept constant at V = 9.4^ = 831 
and we use periodic boundary conditions. According to 
previous simulations [53i (75!, iZfl] the kinetic glass transi- 
tion is around T « 0.435. We analyze here simulations 
at T = 0.15/0.2/0.25/0.30/0.35/0.38/0.40/0.41/0.42 and 
0.43 as they have been described in |4j|. We start 
with 10 independent well equilibrated configurations at 
T = 0.446. Each of these configurations undergoes 
the following sequence of simulation runs. After an in- 
stantaneous quench to T = 0.15 we first run a (NVT) 
simulation for 10^ MD steps. The temperature is kept 
constant by replacing the velocities of all particles by 
new velocities drawn from the corresponding Boltzmann 
distribution every 50 time steps. We then run the sim- 
ulation for 10^ MD steps without the temperature bath 
(NVE) and subsequently raise the temperature instanta- 
neously to the next higher temperature T = 0.2. (NVT) 
and (NVE) runs of 10^ MD steps each follow and the 
temperature is again raised, and so forth. The final con- 
figurations of these (NVT) runs are the initial configura- 
tions for the (NVE) production runs (for 5-10^ MD steps) 
presented in this paper. During each production run the 
positions of all particles (configurations) are stored every 
2000 MD steps which are then used for analysis. 

For the here investigated temperatures the relaxation 
times T are much larger than the waiting time before 
the production runs (at T = 0.446 is r « 8 • 10^). We 
therefore study relaxation processes out of equilibrium 
and find aging effects, which is consistent with previous 



detailed studies of aging of the same binary Lennard- 
Jones system [s^, 01 • 

III. DEFINITION OF JUMP AND JUMP-TYPE 

We focus in this paper on the process of a particle 
escaping its cage, using single particle trajectories ri(t) 
given by the periodically stored configurations. 




FIG. 1: Sketch of a particle trajectory to illustrate the def- 
inition of a jump by comparison of the fluctuations of the 
particle (Ji.ost and the difference in average positions Ar7. 

To distinguish vibrations around an average position 
from a change in the average position, we average 20 
consecutive positions (t) to obtain f7(m) as sketched in 
Fig.n We identify jumps by comparing changes in these 
averaged positions | Af7| with the fluctuations in position 
Ci,ost for each particle i where 

|Ar-| = |r-(m)~f-(r7^-4*20)| (6) 

(see Fig. ^ and (Ti.cst is defined in [t^- We use in spe- 
cific for lAr^l not consecutive but instead average posi- 
tions which are four averages (each of 20 configurations) 
apart. This choice of time interval in |Ar7| is to allow 
identification of not only sudden but also more gradual 
jumps. We define that a jump of particle i occurs when- 
ever 

|Af-p > 20al,, (7) 

Notice that we use a relative criterion, namely for each 
particle i a comparison of |Af7| with its (li^est- Our mo- 
tivation for this relative criterion is that we would like 
to identify jumps of both A and the smaller B particles. 
Also even for particles of the same type their jump size 
might differ due to different cage sizes, crf^^^ is an esti- 
mate for the cage size of each individual particle and is 
therefore used as criterion for the occurrence of a jump. 
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irreversible jump 



reversible jump 



case 1 




case 2 
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FIG. 2: Sketch of typical particle trajectories to illustrate 
the distinction between irreversible and reversible jumps 



IV. NUMBER OF JUMPING PARTICLES 

We apply now the above definitions to identify all 
jumps occurring in our simulations. The number of iden- 
tified jumps depends on the time interval of the produc- 
tion runs and, since the system is out of equilibrium, also 
on the waiting time before the production runs. The 
straightforward counting of jump events would give re- 
versible jumpers (which sometimes jump many times) a 
larger weight. We therefore present the number of jump- 
ing particles. Fig. |2| shows the number of jumping par- 
ticles normalized by the number Na of particles in the 
system of type a £ {A,B}. With increasing tempera- 
ture T the number of jumping particles increases consis- 
tent with an increasing number of relaxation processes. 
A similar temperature dependence has been found indi- 
rectly via the participation ratio in the work 0, \^ |^ . 
More surprisingly, we find that not only the small B- 
particles are jumping but also a considerable fraction of 
A particles. However, the fraction of jumping B particles 
is larger than the fraction of jumping A particles due to 
the smaller size and therefore higher mobility of the B 
particles. Furthermore, both irreversible and reversible 
jumps are not only occurring in a certain temperature 



When we apply the above jump definition, we find 
two types of jumps which we call "irreversible" and "re- 
versible" jumps (see Fig.|2I). A particle which undergoes 
an irreversible jump succeeds in escaping its cage (for the 
time window of the simulation) whereas a particle under- 
going a reversible jump returns back to one of its previous 
average positions. Similar results have been found in pre- 
vious simulations of other systems [1I3. I17 I I21 I [73| . 
However, the present work differs from these that we an- 
alyze all following quantities for the reversible and ir- 
reversible jumps separately. As sketched in Fig. |21 we 
call a jump reversible whenever case 1 or case 2 occurs. 
Case 1 corresponds to the situation of a particle under- 
going multiple jumps and returning to one or more of 
any previous average positions (for details see Sec. Iv)l . 
If two average positions are equal then all jumps which 
happen between the previous position and the recurring 
position are called reversible jumps. In the example of 
Fig. 121 this means in case 1 that all shown jumps are re- 
versible jumps. To increase our resolution in time we use 
not only the information of the averaged positions vlirn) 
but also the complete information of [t] . In case 2 (see 
Fig. 12) the spikes in [t) indicate returns to the average 
position before the jump [80j |. If a jump satisfies neither 
case 1 nor case 2, then it is called "irreversible jump." 




FIG. 3: Number of jumping particles normalized by the num- 
ber of corresponding particle type as a function of tempera- 
ture T. Irreversible and reversible jumpers are distinguished. 



In the following we distinguish between the jump 
events, as they have been defined so far, and the corre- 
sponding jumping particles which may undergo multiple 
jumps of different types. Any reversible jump designates 
the corresponding particle as a reversible jumper for the 
entire time window. 
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FIG. 4: Number of irreversible jumping particles divided 
by the number of both reversible and irreversible jumping 
particles. The inset shows the number of irreversible jump 
events divided by the number of jump events. 



reversible jump irreversible jump 




FIG. 5: This picture illustrates our interpretation of irre- 
versible and reversible jumps and Fig. |1] 

Fig. 01 illustrates the fraction of these jumping parti- 
cles that are irreversible jumpers. At low to interme- 
diate temperatures this ratio is, within the large error 
bars, roughly constant. At intermediate to larger tem- 
peratures irreversible jumpers become more likely than 
reversible jumpers with increasing temperature. This in- 
crease is even more pronounced for all temperatures if 
one considers the number of irreversible jump events (in- 
stead of jumping particles) divided by the number of all 
jump events (see inset of Fig. Gaukel et al. 0,11,01 
come to a similar conclusion via a model for their simu- 
lation data. Their irreversible jumps become more likely 
with increasing temperature. We interpret this increase 
in the fraction of irreversible jumps as sketched in Fig.|Sl 
Both irreversible and reversible jumps start out the same, 
with the jumping particle leaving the cage (formed by the 
neighboring particles) possibly through an opening in the 
cage. In the case of the irreversible jump, the entrance 
of the cage gets meanwhile blocked by a particle, loosely 
speaking the door gets closed, and the jumping parti- 
cle can no longer return and has successfully escaped 
its cage. With increasing temperature all particles be- 
come more mobile, which increases the likelihood of the 
blockage of the entrance back into the cage which in turn 



leads to an increase in the fraction of irreversible jumps 
as shown in Fig. 0] We do not intend to make here a 
statement about the exact geometric process, for exam- 
ple that a door made up of a single particle gets open and 
closed, but more generally the process of rearrangement 
of the cage. Interestingly enough, this blockage happens 
more often in the case of the larger A than the smaller 
B-particles which leads to larger ratios for A than B par- 
ticles in Fig. 01 



V. NUMBER OF AVERAGE POSITIONS 

With each jump a particle either returns to one of its 
former average positions within a cage or to a new over- 
all average position. We call the average positions before 
and after a jump (f7)i and (ri)f (for details about the 
time averages (•)i,t see Fig.^Jand Sec. IVII|I . We use as 
criterion for two average positions to be the same that 
the distance between them A (77) = |(r7)f — (f7)i| and the 
average fluctuations before the jump {af)i (for the defi- 
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number of average positions 

FIG. 6: The number of average positions (as defined in 
Sec. [VJ visited by jumping particles in a) as a function of 
temperature and in b) its distribution for irreversible jump- 
ing B particles. 
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shows an average of the number of distinct average po- 
sitions which are visited by a particle. According to the 
higher mobihty of the smaller B particles, they visit more 
average positions than the A particles. For low temper- 
atures most jumping particles visit only two average po- 
sitions during our simulation. For intermediate to larger 
temperatures, however, not only increasingly more par- 
ticles jump (see Sec. IIV|) but these particles also jump 
more often. The distribution of the number of different 
visited average positions at the highest temperatures (as 
shown in Fig. |^), broadens with increasing T such that 
some particles visit up to seven different average posi- 
tions during the simulation run. 

VI. TIMES 



0-9 A/^ irrev. A 
4—^ A/|^ irrev. B 
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FIG. 7: Sketch to illustrate the definitions of the starting 
time ti and ending time of a jump (here for jumps I and 
II) and the times before the first jump of a particle At\^^^^, 
during a jump Atj, between successive jumps At]-, and after 
the last jump Attail- 

In this section we investigate the time scale of jumps. 
As sketched in Fig. [7|we determine the time duration of 
a jump Aid = if — ii, the time before the first jump of a 
particle Aihcad, the time after the last jump of a particle 
Aitaii, and the time between two successive jumps I and 
II to be Atb = if ^ — 1{ where U and if indicate the start- 
ing and ending time of a jump js^. Notice that Aib is 
only defined if a jump particle jumps twice or more and 
that we analyze Aihcad and A^taii separately as presented 
shortly. For the distinction of irreversible and reversible 
jumps, we assign the jump type of Aib according to the 
jump ending A^b, for example in Fig. [71 the jump type 
of Atb is determined by jump II. This means that A^b 
is a measure of the waiting time before a jump occurs. 
Fig. |S1 illustrates that Aid <C A^b which tells us in hind- 
sight why we could identify jumps as rare events with the 
above described procedure. Since the time resolution is 



FIG. 8: Times between jumps Afb for A (circles) and B 
(diamonds) particles and times during jumps Atd for both A 
and B particles. 



of the order of 1000 (« 20 • 2000 • 0.02), i.e. of the same 
order as Aid, we do not draw further conclusions about 
the temperature dependence of Aid. However, Atb is 
well above our time resolution and small enough to be 
detected during our simulation run of length 1 • 10^. In 
accordance with the picture of reversible jumpers which 
try but do not succeed to escape their cage, these trials 
happen on a shorter time scale than irreversible jumps. 
However, to our surprise, Aib seems to be independent 
not only of the particle type but also of temperature. The 
question arises if this temperature independence of Aib 
is a consequence of the finite time window of our simula- 
tion. To take this finite time interval itot into account, we 
make a crude approximation to correct the time intervals 
Aib, Atjjoad and Attaii by assuming that jumps happen 
equally likely at any time in the window [0, itot] (which 
is not accurate due to aging). The probability Pcorr(Ai) 
of finding the complete interval Ai reduces to the prob- 
ability of finding A< during [0, itot] in the simulation 
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where c is a normalization constant. We may approxi- 



where (-^sim indicates an aver- 



mate c with / ^ **-At 
age over jump events of the simulation. Wc thus obtain 
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Similarly the average times (At)sim may be approxi- 
mately corrected by 
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FIG. 9: Times At\^, Afhoad a^nd Att^il a) for irreversible 
jumps and b) for reversible jumps of A and B particles using 
Eq. dinil. 



Fig-Et and Fig.|3) show the resulting corrected times 
(At)corr for irreversible and reversible jumps respectively. 
For simplification of notation we drop (•)corr in the fol- 
lowing. The tail and head times, A^taii and A^hcad, re- 
veal that we find aging, since A^taii > Athead, i-e. jumps 
are more likely to occur at the beginning of the simu- 
lation than later. This particular aging effect decreases 
with increasing temperature. A^taii and Athcad are both 
larger than A^b because Aib includes only times of par- 
ticles which jump multiple times whereas A^head and 
A^taii include also times of particles which jump only 
once. Interestingly, A^b is (although approximately cor- 
rected) still independent of temperature. Together with 
our results of Secs. lIIII -IVlwe therefore find that with in- 
creasing temperature relaxation processes are accelerated 
via more jumping particles and more multiple jumps, but 
the times between multiple jumps do not become shorter. 
This is contrary to previous results _19j 34, 58j 59, 82] in 
which waiting times decrease with increasing tempera- 
ture. In our system, however, the temperature indepen- 
dence seems to be true not only for the averages but even 
for the distribution Pcori (Aib) (see Fig. I1U|I . We inter- 
pret this temperature independence as being due to aging 




FIG. 10: Log-log plot of P(Afb) using Eq. @ for irreversible 
jumps of A and B particles. The line is a linear fit to the data 
with slope —0.84. 



which is consistent with our data for Att a.ii and Athoad 
and with the results of Doliwa et al. [53, |53 who find 
that the distribution of waiting times initially is temper- 
ature independent and becomes temperature dependent 
at later times. As shown in Fig. ^| for irreversible jumps 
of A and B particles, P(Aib) approximately follows a 
power law P{Ath) oc At^'^ with v = 0.84. Even though 
we expect the system to show normal diffusion on a time 
scale much larger than our simulation run, the intermedi- 
ate dynamics is subdiffusive ^8^. Please notice however 
that a direct comparison of our result with the scheme 
presented in Fig. 1 of [s^l is not possible because in our 
case the multiple jumps of one particle at different times 
as well as the jumps of different particles are not inde- 
pendent, which changes the dynamics ::84i. 
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FIG. 11: Sketch for the definitions of A_Rif, Ai?max and 

Aiiavg. 



VII. JUMP LENGTHS 

In the last sections we have investigated how many 
particles jump and how often jumps occur. Next we 
study how far these particles jump. To be able to test if 
our qualitative results are dependent on the definition of 
jump length, we use three quantities: Ai?if, Ai?rnax and 
Ai?avg as sketched in Fig. ^2 For a jump of particle i 
the jump distance of the jump startin g at ti and ending 
at tf (for the definition of tj and if see [S^) is 

AiJif- |rr(ti)-rr(if)| (11) 

and the maximal distance being detected during the jump 

Ai?max = max |fr(t) - vl{t - 3200) | (12) 

(here 3200 = 4 • 20 • 2000 • 0.02 as in ^ and Sec. ^ 
where t is varied over times t which satisfy ii,detect < ^ < 
if,detect (US)- The third length Ai?avg is less dependent 
on fluctuations than Ai?if and Ai?max. It is the distance 
of the overall average positions before and after the jump 

Ai?avg = |(fr)f-(fr>il (13) 

For the averages (•)t and the positions Vi{t) are av- 
eraged over times before and after the jump excluding a 
broadened time window around the jump [t{ — 800, if + 
800]. Fig. 112b shows that the different jump lengths have 
the same qualitative behavior, namely, as one might ex- 
pect, an increase with increasing temperature T. We 
therefore show in the following mainly results for only 
Ai?avg but find very similar behavior in the case of Ai?if 
and Ai?max- An increase in jump size with increas- 
ing temperature has been seen in previous simulations 
[1 E iHl H Ip. However in the case of 0, H 
the jump size is for a jump in the MSD, i.e. an average 
over all particles, which includes two effects: an increas- 
ing number of jumping particles and an increase in the 
jump size of single particle jumps. Our approach has the 
advantage of keeping these two effects separate. Further- 
more, in our case the single particle jumps are larger than 
in HIHIilES 

Fig.[T% is the same as Fig.lT^ for Ai?avg but broken 
down by both jump and particle type. The smaller B par- 
ticles are jumping further than the A particles. Further- 
more reversible jumps are shorter than the irreversible 
jumps, which is due to two features in the distribution 

^"(Ai?avg). 

As illustrated in Fig.Elboth irreversible and reversible 
jumps have a peak at about 0.8 and 1.0 for A and B 
particles respectively. These peaks are relatively similar 
for reversible and irreversible jumps, taking aside that for 
irreversible jumps this peak position is slightly shifted to 
the right and slightly broadened (partly due to multiple 
jumps which are not resolved in time). This similarity 
and the peak positions around unity are consistent with 
our picture that irreversible and reversible jumps start 
out the same, namely with a particle jumping out of its 
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FIG. 12: a) The jump sizes AfJjf, A_Rmax and A_Ravg 
(sketched in Fig. 1111 and defined in the text) as function of 
temperature T. The averages are over both irreversible and 
reversible jumps of both A and B particles, b) The jump size 
A7?avg as function of temperature T separately for irreversible 
and reversible jumps of A and B particles. 

cage. A new feature, however, is that the distribution 
of reversible jumps for T — 0.43 (and similarly for all 
other temperatures) has a bimodal distribution with an 
additional peak at a distance smaller than 0.2, contrary 
to the distribution of irreversible jumps. This seems to 
indicate for reversible jumps another process than a jump 
out of a cage because the latter would be of the order 0.5 
and larger. 

The increase of Ai?avg with increasing temperature is 
a consequence of both a shift of the peak position and a 
broadening of the distribution of jump sizes. This is illus- 
trated in the inset of Fig. ^| for irreversible jumps of A 
particles and we find similar distributions for irreversible 
jumps of B particles. At the largest investigated tem- 
perature T = 0.43 some of the particles move as far as 
four particle spacings. Some of these large jumps might 
correspond to multiple jumps of a smaller time window 
than our time resolution. 
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FIG. 13: The distribution of jump sizes A7?avg for irreversible 
and reversible jumps of A and B particles at T = 0.43. The 
inset shows Af?avg of irreversible jumping A particles for dif- 
ferent temperatures. 
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FIG. 15: Fluctuations in position a for jumping particles 
Cjumpj of the original estimate for jumping particles (Test, and 
for non-jumping particles (Jnojump- For the average of ajump 
fluctuations during the jump have been excluded from the 
average. 




FIG. 14: Log-log plot of the distribution of jump sizes A/2avg 
for irreversible jumps of A and B particles including all tem- 
peratures. The lines are linear fits with slopes —6.3 for A and 
—4.8 for B particles. 

For irreversible jumps of size unity and larger we find 
that the distribution for all temperatures follow roughly a 
power law P{AR) « AR'^ (see Fig.[nil with ^ = 6.3 for 
A and fi — 4.8 for B particles. This is in accordance with 
subdiffusive behavior for intermediate times as presented 
in Sec. I VII for the distribution of waiting times js^]. 

As commented on earlier (see Sec. U and Sec our 
definition for a jump is not based on a specified size but 
instead a multiple of the fluctuations of the particle. Let 
us therefore next look at the fluctuations. Fig.^lshows a 
comparison of the fluctuations of jumping particles CTjump 
and of non-jumping particles Unojump- Also included in 
Fig-Elis an average over all particles (Test of the estimates 
for fluctuations Ci^est as they have been used for the jump 



identification (see Sec. IIIII and [t^). (Test is very similar 
to 

fnojump and represents the fluctuations of an average 
particle. For the average of crjump we exclude fluctuations 
during the jump: for jumping particle i we take 



1/2 



' i,jump 



(14) 



where cr^ and (•)i_t are defined as in |73 and earlier in this 
section respectively. Although jump times are excluded, 
the fluctuations of jumping particles are clearly larger 
than the fluctuations of non-jumping particles. This 
means that jumping particles are not only moving far- 
ther than an average particle during single events but 
are also oscillating within their cages with a larger than 
average amplitude. 

Fig. \Wh shows an average of jump size relative to the 
particle's fluctuation Ai?i_avg/o'ijump for jumping parti- 
cles i. Ai^avg/cjump seems independent of temperature 
for reversible jumps and for irreversible jumps slightly 
increasing with increasing temperature. 
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AR /a , 

avg est 

FIG. 16: a) Normalized jump size AT^avg/ c^jump ^is function of 
temperature separately for irreversible and reversible jumps 
of A and B particles, b) Distribution of normalized jump 
size P(A_Ravg/o"ost) of irreversible jumping A particles and 
reversible jumping B particles in the inset. 

The distribution of AR/a gives us the opportunity to 
estimate the influence of the cutoff in our definition of 
a jump. Similar to our approach for the search of a 
jump, we use (Test) and obtain P(Ai?avg/o'est) as shown 
in Fig. 116b for irreversible jumps of A particles. We find 
a peak at (AiJavg/cest) ~ 5.5 (for B-particles ~ 6.5), 
i.e. larger than our cutoff at (AR/a) = y/20 « 4.5. This 
gives us the hope that we have included the major contri- 
bution of jumps consistent with our approach. We have 
to bear in mind, however, that Ai?avg/o'cst is not the 
identical measure to the one used in our search proce- 
dure. Furthermore, in the case of reversible jumps there 
is some additional peak at smaller values of (Ai?avg/o'cst) 
(see inset of Fig.lTBb) which is consistent with Fig.ll3land 
might indicate some additional process. 



VIII. ENERGY 

In the previous sections we have investigated relaxation 
processes by solely using the positions (trajectories) of 



single particles. As the work [ssllsllSTllsllsglleot 

shows, an approach via potential energies can also be 
fruitful. We therefore study in this section what happens 
to the potential energy during jumps. As sketched in 




FIG. 17: Sketch for the definitions of {E)i and {E){. 



Fig. El we still use single particle trajectories to identify 
the jumping particle and do similar averaging for the 
according time windows and (•)t (see Sec. IVII|I . We 
investigate both the total potential energy per particle 

1 ^ 

i=l j>i 

and the single particle potential energy of particle i 

E,{t) = Y^Vo,p{n,{t)) (16) 

where VQ^(ry (t)) is defined in Eq. ©. To obtain in 
addition the energies of the inherent structures [s^ we 
minimized the configurations ({ri(i)}) via a steepest de- 
cent procedure. With the thus obtained configuration 
({r°(i)}) we determine 

1 ^ 

and 

Ent)-J2^Mr^jit)) (18) 

Similar to the jump size in position (see Fig. 1171 and 
Sec. IVIIjl we then determine 

ASavg = {Eit))f - {E{t))i (19) 

and 

AEu = E{U) ~ E{U) (20) 
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and similarly for the energies of Eq. (|16|l - H18|) . 

We find in most of our analysis that AiSavg and AE^if 
(and all other equivalents for E^.Ei and E^) are show- 
ing the same qualitative behavior, and therefore present 
in the following results for Ai?avg using Eq. p9|l and 
□ □ 
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FIG. 18: As a function of temperature T a) jump size in total 
potential energy Ai?avg and b) absolute value of the jump size 
in minimized single particle potential energy |Ai5favgl- 

Fig. 118b shows Ai^avg as a function of temperature av- 
eraged separately over irreversible and reversible jumps 
of A and B particles. Notice that the energy jumps are 
small compared to energy values {E)-^^ « —7. The ir- 
reversible jumps lower the energies more than reversible 
jumps and more so for A than B particles, because a 
jump of a larger A particle results in a bigger change in 
environment than for a smaller B-particle. The lowering 
of the total potential energy happens mostly at inter- 
mediate temperatures and less at the lowest and highest 
investigated temperatures. This might be an aging ef- 
fect: at low temperatures the system is basically frozen 
in, at intermediate temperatures the system has time to 
partially age, and at higher temperatures the system has 
already aged before the production run starts. This ef- 
fect is lost when we look at the jumps in single particle 



energies AE'i.avgC?^) which is basically zero for all jumps. 
We therefore find that the total potential energies Ai?avg, 
AiJif, Ai?avg, and Ai?j° show a more systematic depen- 
dence on temperature than their single particle equiva- 
lents ASj^avg etc.. For the absolute values |Ai?avg| etc. 
and all following quantities, however, we find the oppo- 
site, i.e. the single particle equivalents |Ai?i^avg| etc. 
show a more pronounced behavior. We therefore show 
from now on the results which use Eq. H16|l and Eq. H18() 
only. 

We next look at the absolute value |Ai??,j^g| as a func- 
tion of temperature (see Fig. 118b). As one might ex- 
pect, the irreversible jumps show larger changes in en- 
ergy than reversible jumps, due to larger jumps in posi- 
tion fFig. ll2b ) and therefore more change in the environ- 
ment of the jumping particle. And similarly the larger 
A particles in comparison with the smaller B particles 
experience a larger change in energy because A particles 
are surrounded by more neighbors, i.e. a larger environ- 
ment. With increasing temperature the absolute value of 
the jump in energy fFig. ll8b ) is increasing (with the one 
exception of irreversible B particles) which is consistent 
with the increase in jump position fFig. 112( 3). 

We define the fluctuations in energy at time t = 
m ■ 2000 • 0.02 



\ I \ m'— m — 9 / ^ 

+ (yo E" {E.im') -E.irn')A I) 

\ m'—m — 9 ' f J / 



-1/2 

(21) 



similar to the fluctuations in position |79l | where Ei[m) is 
m+10 

the time average ^ ^ Ei{m'). The energy fluctua- 

m'— m — 9 

tions arc also increasing with increasing temperature and 
are larger for A than B particles (see Fig. 119b ) . The fluc- 
tuations of irreversible and reversible jumps are however 
very similar, which is consistent with the picture that re- 
versible jumps are "failed" irreversible jumps and in that 
sense start out the same way as irreversible jumps. 

Fig. 119b is a comparison of energy fluctuations of jump- 
ing (both irreversible and reversible and of both A and B) 
particles ijR . and of non- jumping particles ctr 
|87l |. In o'E.^^^p we exclude the times during the jump 
by using the time windows We find that the en- 

ergy fluctuations of jumping particles are larger than for 
non-jumping particles which is consistent with the corre- 
sponding fluctuations in position (see Fig. 115(1 . 



By normalizing jAi^i avgl with the fluctuations (jEi we 
obtain Fig. 120b .. Similar to the case of normalized jumps 
in position (Fig. Et) we find that | Ai?i_avg|/cBi is for 
reversible jumps basically independent of temperature. 
For irreversible jumps, however | Ai^i^avgl/cfii decreases 
significantly with increasing temperature. This indicates 
that with increasing temperature the jumps are increas- 
ingly more driven by fluctuations. Contrary to the equiv- 
alent in position (Fig. 116b ) , the normalized energies are 
signiflcantly smaller, in the range 0.5 - 1.4 rather than 3 
- 7, and the distribution of lAiJi^avgl/fEi (see Fig.lJOt) is 
monotonous decreasing with a less far reaching tail than 
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FIG. 19: Fluctuations in single particle potential energy as a 
function of temperature. The average is in a) over irreversible 
and reversible jump events of A and B particles separately and 
in b) over jumping particles {as. jump) ^-^id over non-jumping 
particles {as. ■ )■ 



IX. CONCLUSIONS 

We study the dynamics of a binary Lennard-Jones sys- 
tem below the glass transition. Our focus is on jump- 




\AE.\/cy^ 



FIG. 20: Absolute value of single particle potential energy 
divided by its fluctuation \ AEi^avg\/o'Ei in a) as a function of 
temperature and in b) its distribution . 



processes, which we identify via single particle trajecto- 
ries. Two kind of jumps are found: "reversible jumps," 
where a particle jumps back and forth between one or 
more average positions, and "irreversible jumps," where 
a particle does not return to any of its former average 
positions, i.e. successfully escapes its cage of neighbors. 
Both irreversible and reversible jumps of A and B parti- 
cles occur at all temperatures. With increasing temper- 
ature more particles jump, more average positions are 
visited, the jump size both in position and in the ab- 
solute value of the potential energies \AE\, \AE°\, \AEi\ 
and |Ai?°| (total and single-particle, not minimized and 
minimized) increase, and the fluctuations in position an 
and potential energy as increase. The fluctuations are 
larger for jumping particles than for non-jumping par- 
ticles even if fluctuations during the jump are not part 
of the average, which indicates that jumping particles 
are not only during jump times more mobile than non- 
jumping particles. The ratio \AEi\/aEi of irreversible 
jumps decreases with increasing temperature. This con- 
firms a commonly used assumption that with increasing 
temperature the irreversible jumps become more driven 
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by fluctuations. 

With increasing temperature also irreversible jumps 
become proportionally more frequent than reversible 
jumps. We interpret this such that irreversible and re- 
versible jumps are similar in that sense that a particle 
tries to escape its cage. In the case of reversible jumps 
the particle finds its way back into the cage whereas in 
the case of an irreversible jump the path back into the 
cage becomes blocked due to rearrangements of the cage. 
At larger temperature these rearrangements of the cage 
become more likely (since for example the fluctuations in- 
crease) and therefore irreversible jumps occur more often. 
Irreversible and reversible jumps show in most quantities 
qualitatively the same behavior, such as their temper- 
ature dependence of jump size in position and energy, 
and differ only in size. An exception to this similarity is 
the distribution of jump sizes, which suggests that there 
might be an additional jump process in the case of re- 
versible jumps indicated by an additional peak at small 
jump sizes. 

The most surprising result of our work is that the times 
between successive jumps are independent of tempera- 
ture. This is most likely due to aging, which could mean 
on the time scale of our simulation that A<b reflects the 
time scale of T = 0.446 from which we quenched the sys- 
tem, and that Atb might show temperature dependence 
at later times. The latter would be consistent with the 
work of Doliwa et al. [sM l59j| who find in their simu- 
lations of very long times that the distribution of wait- 
ing times (however, identified via the collective quantity 
of the minimized potential energy) initially is tempera- 
ture independent and of different power law than at later 
times when temperature dependence occurs. We do find 
aging in that sense that the times before the first jump 
are shorter than the times after the last jump of a particle 
(an effect which becomes less for temperatures near Tc). 



Future work (which is in progress) will tell us if in our 
system Aib becomes temperature dependent after longer 
times. However, the dynamics does depend on tempera- 
ture for all other here investigated quantities such as the 
number of jumps and the fraction of reversible jumps. 

Another question which wc would like to raise for fu- 
ture work is how much single particle jumps can tell us 
about relaxations of a glass which might include collec- 
tive motion. It seems plausible that single particle jumps 
are strongly correlated with collective jumps and in that 
sense give us very similar information to studies of collec- 
tive quantities. We find that many of the single particle 
jumps in our system are spatially and temporally cor- 
related, i.e. are showing collective motion. Systematic 
studies for confirmation remain to be done in the future. 
We also find the signature of many-particle effects con- 
firmed in the result of the increasing ratio of irreversible 
jumps with increasing temperature, since this indicates 
that not only the jumping particle itself but also its cage 
are dependent on temperature. There might, however, 
be also collective jumps, where each particle jumps a too 
small amount to be detected by our search algorithm. 
It remains to be studied of how much importance such 
processes are in our system. 
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